ACCEPTED BY: The Astronomical Journal 

Preprint typeset using I^-T^X style cmulatcapj v. 08/13/06 



O 

(N 

X> 
IX, 



O 

u 

43 
6 

03 



HIERARCHICAL STRUCTURE FORMATION AND MODES OF STAR FORMATION IN HICKSON COMPACT 

GROUP 31* 

S. C. Gallagher, 1 ' 2 P. R. Durrell, 3 D. M. Elmegreen, 4 R. Chandar, 5 J. English, 6 J. C. Charlton," C. 



Gronwall, J. Young, P. Tzanavaris 



C. 



Mendes de Oliveira, 
Iz Ann Zabludoff 14 



B. Whitmore, 5 A. E. 



8,9 K. E. Johnson 

HORNSCHEMEIER, 8 APARNA MAYBHATE 13 

accepted by: The Astronomical Journal 

ABSTRACT 

The handful of low-mass, late-type galaxies that comprise Hickson Compact Group 31 is in the 
midst of complex, ongoing gravitational interactions, evocative of the process of hierarchical structure 
formation at higher redshifts. With sensitive, multicolor Hubble Space Telescope imaging, we charac- 
terize the large population of < 10 Myr old star clusters that suffuse the system. From the colors and 
luminosities of the young star clusters, we find that the galaxies in HCG 31 follow the same universal 
scaling relations as actively star-forming galaxies in the local Universe despite the unusual compact 
group environment. Furthermore, the specific frequency of the globular cluster system is consistent 
with the low end of galaxies of comparable masses locally. This, combined with the large mass of 
neutral hydrogen and tight constraints on the amount of intragroup light, indicate that the group is 
undergoing its first epoch of interaction-induced star formation. In both the main galaxies and the 
tidal-dwarf candidate, F, stellar complexes, which are sensitive to the magnitude of disk turbulence, 
have both sizes and masses more characteristic of z = 1-2 galaxies. After subtracting the light from 
compact sources, we find no evidence for an underlying old stellar population in F - it appears to be 
a truly new structure. The low velocity dispersion of the system components, available reservoir of 
H I, and current star formation rate of ~10 M Q yr _1 , indicate that HCG 31 is likely to both exhaust 
its cold gas supply and merge within ~ 1 Gyr. We conclude that the end product will be an isolated, 
X-ray-faint, low-mass elliptical. 

Subject headings: galaxies: star clusters — galaxies: interactions — galaxies: evolution — galaxies: 
clusters: individual (HCG 31) 
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1. INTRODUCTION 

Massive elliptical galaxies are preferentially found in 
high density environments such as galaxy clusters - 
galaxy morphology and environment are clearly cou- 
pled in the local Universe (Dressier 1980). However, 
the mechanism for morphological transformation from 
star-forming and disk-dominated to quiescent and bulge- 
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dominated is still uncertain. For example, is ram- 
pressure stripping important as a galaxy plows through 
the hot intracluster medium (ICM)? Alternatively, grav- 
itational interactions from harassment to mergers might 
dominate the transformation process. Secular evolution 
whereby gas is sheparded inwards, thus triggering star 
formation and exhausting the material for subsequent 
generations of stars, could be important. The action of 
an active galactic nucleus (AGN) may also significantly 
alter its host in dense regions. 

Furthermore, the question of when morphological 
transformation occurs in the hierarchical build-up of a 
galaxy cluster from subunits of galaxy groups is still 
open. Already, the demographics of cluster galaxies as 
a function of redshift indicate that some galaxy evolu- 
tion occurs prior to cluster infall (e.g., Dressier et al. 
1997). Determining if galaxies are primarily processed 
while still contained within the unit of the group will shed 
light on the influence of the ICM (e.g., Zabludoff et al. 
1996; Zabludoff & Mulchaey 1998, 2000). Studies of com- 
pact group galaxies offer the promise to further elucidate 
this issue, as in their cores they have the same galaxy 
number densities as galaxy clusters, but without (with 
a few exceptions) an X-ray emitting intragroup medium 
(ICM). The lower velocity dispersions (a ~ 10 2 kms -1 ) 
in groups compared to clusters (a ~ 10 3 kms^ 1 ) also 
prolong gravitational interactions. This context is rel- 
evant therefore for evaluating the role of gravitational 
interactions in morphological transformation before the 
cluster potential or ICM becomes dominant. 

Recent Spitzer results on compact group galaxies sug- 
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gest that the transformation from neutral gas-rich and 
actively star-forming galaxies to neutral gas-poor and 
quiescent galaxies can occur quite rapidly in this envi- 
ronment. Specifically, Johnson ct al. (2007) found that 
most member galaxies in compact groups in an early evo- 
lutionary state (defined by a high ratio of their H I to 
dynamical masses) are actively star-forming as judged 
from their infrared luminosities and colors, while those 
with little H I are dominated by the stellar photosphcric 
emission of old stellar populations. While this is perhaps 
not surprising, there is a notable gap in the distribution 
of infrared spectral indices of the compact group popu- 
lation compared to a sample of local galaxies matched in 
optical luminosity from the SINGS sample (Kennicutt et 
al. 2003); the SINGS galaxies have a continuous distribu- 
tion of infrared properties (Gallagher et al. 2008). The 
lack of compact group galaxies in the intermediate re- 
gion of infrared color-space suggests that this location is 
crossed quickly (Johnson et al. 2007; Walker et al. 2009). 
This picture of galaxy processing in this environment is 
further demonstrated by the H I deficiency of late-type 
compact group galaxies relative to those in loose groups 
and the field (Williams & Rood 1987; Verdes-Montenegro 
et al. 2001). 

Local Hickson Compact Groups (HCGs; Hickson 1993) 
are therefore promising testbeds for investigating mor- 
phological galaxy transformations nearby that are di- 
rectly relevant to galaxy evolution from 'blue cloud' 
to 'red sequence' in the color-absolute magnitude plane 
(e.g., Hogg et al. 2004; Balogh et al. 2004; Willmer et al. 

2006) . In currently dense environments (i.e, galaxy clus- 
ters), this change likely occurred at z = 1-2 (e.g., Cuc- 
ciati et al. 2006; Franzetti et al. 2007; Cooper et al. 

2007) . Of the nearest HCGs, HCG 31 (NGC 1741) is 
perhaps the best example of catching this process in 
progress. This system is comprised entirely of low mass 
(1-8 x 10 9 M Q ; Verdes-Montenegro et al. 2005), late- 
type galaxies in a compact (~ 1.4' [~ 23 kpc] diameter) 
configuration. It is clear from the multicolor image in 
Figure 1 that each of galaxies in the group is disturbed, 
and tidal tails, bridges, and dwarfs are prominent. 

In this paper, we present deep, three band HST ACS 
images of HCG 31, using high angular resolution and 
sensitivity to probe modes of star formation from the star 
cluster (~pc) to star complex and dwarf galaxy (~kpc) 
scale. Our goal is to investigate hierarchical structure 
formation through an in-depth study of this system. In 
particular, we want to study the mechanism by which 
neutral gas is consumed in dense and dynamically young 
environments. 

1.1. Studies of HCG 31 to Date 

From visual inspection of Figure 1, the primary galax- 
ies that comprise HCG 31 (A+C, B, and G) are clearly 
disturbed, and the presence of bright tidal structures (E, 
F, and H) further supports the importance of gravita- 
tional interactions in the recent history of the system. 
An additional galaxy, Q to the north, is shown to be 
part of the group based on radial velocity (Rubin et al. 
1990) and H I structure (Verdes-Montenegro et al. 2005) 
that link it to other group members. UV, optical, and 
infrared photometry of A+C indicate active recent and 
ongoing star formation (e.g., Iglesias-Paramo & Vilchez 
1997; Johnson et al. 1999; Johnson & Conti 2000; John- 



son et al. 2007). This conclusion is supported by spectro- 
scopic detections of Wolf-Rayet signatures in A+C and 
F (Kunth & Schild 1986; Rubin et al. 1990; Conti 1991; 
Lopez-Sanchez et al. 2004). 

The high levels of star formation and irregular mor- 
phology of the A+C complex imply a recent starburst 
triggered by their mutual interaction. Further study of 
the H I morphology and velocity field of the entire sys- 
tem indicates that not just one but multiple interactions 
have led to the current configuration (e.g., Amram et al. 
2004; Mendes de Oliveira et al. 2006). Notably, much of 
the neutral hydrogen in the group has been removed from 
the galaxies (Verdes-Montenegro et al. 2001) yet remains 
bound to the group in tidal structures; 60% percent of 
the H I is located in four tidal tails and one bridge link- 
ing six of the member galaxies (Verdes-Montenegro et al. 
2005). In sum, though the group as a whole is not H I 
deficient, the individual galaxies are (Verdes-Montenegro 
et al. 2001). 

While Rubin et al. (1990) conjectured that HCG 31 
will coalesce into a single large elliptical in only a few 
orbital periods, this is still an open issue. Specifically, 
will HCG 31 transition in the near term (within ~ 1 Gyr) 
into an evolved state such as a 'fossil group' - a single, 
elliptical galaxy with a bright X-ray halo (Jones et al. 

2003) ? As HCG 31 is one of the most compact groups 
in the Hickson et al. (1992) Atlas with most of the H i 
stripped from its individual galaxies, it is perhaps the 
most likely late-stage merger candidate (Williams et al. 
1991). 

Can the modes of star formation — from the pop- 
ulations of parsec-scale star clusters to kpc-scale star- 
forming complexes and dwarf galaxies — in this compact 
galaxy group elucidate the mechanism of galaxy trans- 
formation in the high redshift universe? With our new 
ACS imaging as well as IR+UV constraints on the cur- 
rent star formation rates (SFRs) throughout the system, 
we aim to shed new light on the current state of star 
formation in HCG 31 as well as its future evolution. 

At the H I redshift, z = 0.01347+0.00002 (Meyer et al. 

2004) , of HCG 31, its distance is 58.3 Mpc (yielding a 
distance modulus, DM =33.83) and 1" corresponds to 
275 pc (CIm — 0.3; Q\ = 0.7; h wo = 0.7; Spergel et al. 
2007). Using a Galactic E(B - V) = 0.051 (Schlegel 
ct al. 1998), all magnitudes and colors throughout the 
text have been corrected for Galactic reddening. 

2. OBSERVATIONS 

Observations of HCG 31 were obtained with HST 
ACS/WFC using the F435W, F606W, and F814W fil- 
ters. Throughout, we use B 4S5 , V^, and I 81 4 to re- 
fer to the F435W, F606W, and F814W ACS filters and 
the magnitudes measured in those filters, respectively. 
The data were taken on August 8, 2006 with expo- 
sure times of 1710 seconds in F435W, 1230 seconds in 
F606W, and 1065 seconds in F814W. The observations 
for each filter were taken with three equal exposures, us- 
ing a three-point dither pattern. The implementation 
of MultiDrizzle 16 in the HST/ ACS data pipeline pro- 
vides combined, geometrically corrected, and cosmic ray 

16 MultiDrizzle is a product of the Space Telescope Sci- 
ence Institute, which is operated by AURA for NASA. Sec 
http://stsdas.stsci.edu/pydrizzle/multidrizzle. 



Modes of Star Formation in HCG 31 



3 



cleaned images. For the analysis of point sources, we 
used the standard HST pipeline products with a nominal 
pixel scale of 0'.'05 per pixel. For analysis of the extended 
sources, we ran MultiDrizzle with the pixel scale set to 
0'.'03 per pixel to improve the spatial resolution. 

In addition to the new HST observations that are the 
primary focus of this paper, we also present new Swift 
and GALEX observations of HCG 31, and incorporate 
the Spitzer IRAC and MIPS images presented in John- 
son et al. (2007). This is part of a long-term, in-depth 
multiwavclcngth investigation into a sample of 12 HCGs 
using the combination of HST three-color imaging (six 
complete), Chandra ACIS imaging spectroscopy (eight 
total as part of archival and new observations), Spitzer 
IRAC and MIPS imaging (already in hand; Johnson et al. 
2007), and Swift Ultraviolet Optical Telescope (UVOT) 
and GALEX multifiltcr imaging (Tzanavaris et al., in 
prep.). 

For the purposes of associating star-forming structures 
within the individual features in HCG 31, we have de- 
fined polygonal boundaries for the galaxies and tidal re- 
gions as marked on the three-color ACS image shown in 
Figure 2. 

2.1. Star Cluster (SC) Candidate Identification and 
Photometry 

At the redshift distance of HCG 31 of 58.3 Mpc, star 
clusters are expected to be unresolved or marginally re- 
solved with a 0'.'05 ACS pixel corresponding to 13.8 pc. 
Therefore, as a first step to finding star cluster candi- 
dates, the entire ACS field was searched for point sources. 

Point-source photometry was performed on the 
pipeline-drizzled images using the point source function 
(PSF) fitting algorithm of Stetson (1987) as implemented 
in IRAF. 17 Given that the Vqqq image was the deepest of 
the three, it was the primary one used for point-source 
identification and filtering. Magnitudes arc given in the 
Vega system using the zeropoints appropriate for August 
2006: 18 26.406, 25.767, and 25.520, for V 606 , S 435 , and 
Isi4, respectively. 

We searched for point-source candidates using the 
median-divided image with a threshold of 0.25 

counts; source- free background regions have standard de- 
viations of ~ 0.07 counts. The weight-map image from 
the pipeline processing was used to filter the candidate 
list to remove artifacts from chip edges and the diffrac- 
tion spikes from saturated stars; point source candidates 
were rejected if their weight map values were less than 
the median minus twice the standard deviation of the dis- 
tribution. Visual inspection confirmed that this method 
effectively screened out spurious detections without dis- 
carding real point sources. 

PSF models in each filter were constructed from the 
brightest, isolated, and unsaturated stars. Their radial 
profiles were checked for irregularities, and 13, 13, and 18 
PSF stars were used in Vqqq, -B435, and Isu, respectively. 
Moffat functions were used to fit the PSF in each case, 
with a constant PSF in Vgoe, and linearly variable PSFs 

17 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 

18 http://www.stsci.edu/hst/acs/analysis/zeropoints/ 



for -B435 and Isi4- The daophot photometry was com- 
pared to aperture photometry for the PSF stars within 
a 0'.'5 (10 pix)-radius circular aperture in order to deter- 
mine the aperture correction. To these values, the aper- 
ture correction from 0'.'5 to infinity from Sirianni et al. 
(2005) was added. Finally, the photometry was corrected 
for Galactic extinctions of Aqqq=Q. 144, ^435=0. 210, and 
A 8 i4=0.093. 

For the final star cluster (SC) candidate catalog, ob- 
jects were required to have My 606 < — 9 (to eliminate con- 
tamination from individual young stars), be detected in 
all three filters with magnitude error <7 mag < 0.3 and 
-2< sharp <2 (sharp is a measure of the width of the 
data relative to the model psf; sharp = indicates a 
width well-matched to the psf) in each and x < 3 in 7 8 i 4 . 
We used Isi4 for the \ criterion because the model PSF is 
best determined in that image:. Furthermore:, if SC candi- 
dates are marginally resolved because of extended, nebu- 
lar Ha emission, this will not be a factor in 7 8 i 4 . (We do 
not want to eliminate these interesting candidates from 
our catalog.) Finally, color cuts of -B 4 35-V606< 1-5 and 
^606 _ -^8i4< 1-0 were imposed to reduce contamination 
from foreground Galactic stars. These cuts were chosen 
to eliminate a distinct clump of sources in color-color 
space that are well-isolated and far from the expected 
colors of star clusters. The SC candidate numbers by re- 
gion are given in Table 1, and a -B435-V1506 vs - ^606 _ ^8i4 
plot with the sources that meet our criteria as outlined 
above are presented in Figure 3. To identify the youngest 
clusters with nebular emission, we draw a diagonal line 
in Figure 3 that crosses the evolutionary tracks with neb- 
ular emission at r ~ 5 x 10 6 yr. This conservative line 
takes into account typical photometric errors such that 
all SC candidates to the left of that line have colors sig- 
nificantly distinct from the evolutionary tracks without 
nebular emission. The numbers and fractions of SC can- 
didates with these colors for each region of HCG 31 are 
indicated in columns 5 and 6 of Table 1 . 

We imposed stricter point-source selection criteria to 
identify globular cluster candidates, and extend the mag- 
nitude limit to Mv 606 < —7.8 to reach the expected peak 
of the globular cluster luminosity function. Following 
Rejkuba et al. (2005), we used a hyperbolic filter in the 
sharpness vs. magnitude plane and a hyperbolic filter 
in the magnitude error vs. magnitude plane set using 
the output photometry from the V606-band complete- 
ness tests (see Fig. 3 in Rejkuba et al. 2005 and §2.2) 
to include ~ 97% of the detected artificial point sources. 
Point sources were also required to have cr mag < 0.3, 
X < 3, and —2 < sharp < 2 in all three filters. 
These criteria rejected a large fraction of the brightest 
(Mv 606 < —9) and bluest (in V606~-^8i4) point sources, in- 
dicating that many of them are marginally resolved. This 
is perhaps not surprising, as the wide Vgo6 filter includes 
Ha emission, and most of these sources have colors con- 
sistent with a strong nebular contribution. The spatial 
scale of the Ha-emitting region surrounding a young star 
cluster is frequently larger than that of the stars alone. 
Those candidates with colors consistent with older star 
clusters are typically fainter and unresolved, as expected. 

2.2. Completeness 

Using addstar in the daophot package, we added 3000 
artificial stars to the pipeline-drizzled Vqo6 image evenly 



4 



Gallagher et al. 



distributed between magnitudes 24 through 28 (before 
aperture and extinction corrections), and then ran the 
photometric pipeline on these data. Sources were consid- 
ered detected if their input and output positions matched 
within 0'.'075 (1.5 pix). Using the distributions of sharp- 
ness, sharp, vs. magnitude, and magnitude error, cr mag , 
vs. magnitude, we then set the same strict point-source 
selection criteria described in §2.1. Of the detected 
sources, these filters eliminated an additional ~ 3%. 
From the ratio of numbers of detected point sources to 
input artificial stars as a function of magnitude, the Veoe 
90% and 50% complete magnitude limits were found to 
be 26.43 and 27.42 (including aperture and extinction 
corrections), respectively. For the B435 and 7 8 i4 tests, 
the same input artificial star catalog was used as for the 
V606 image with the positions corrected for the image off- 
sets. The V<306 detected objects were then run through 
the -B435 and Isu photometric pipelines. Objects were 
filtered using only the constant absolute <7 mag , X; an d 
sharp limits defined above, which eliminated 1 — 2% of 
the sources. The subsequent -B435//814 90% and 50% 
completeness limits were 26.73/26.51 and 27.33/27.18, 
respectively. Given the different aperture and extinction 
corrections in each filter, the input artificial star cata- 
log had effective colors of -E?435-V606= —0.13 and Vqoq- 
J 8 i 4 = 0.04. 

Given that most of the young SC candidates are lo- 
cated within the group galaxies, we also determined the 
completeness limits within A+C and B, with boundaries 
corresponding to the regions shown in Figure 2 where the 
average surface brightness and background noise are con- 
siderably higher. The 90% V606 / B435 / -^814 completeness 
limits are 25.02/24.96/25.10. For a distance modulus of 
33.83, even these limits are significantly fainter than the 
apparent magnitude of 24.83 corresponding to an abso- 
lute V606 luminosity limit of My 606 < —9. 

2.3. Star-Forming Complexes 

Complexes are isolated, high surface-brightness, star- 
forming knots whose sizes and masses are sensitive to 
the velocity dispersion in the gas. Visually, they are of- 
ten identified as "beads on a string" along spiral arms. 
Galaxy disks stirred up by tidal interactions would be ex- 
pected to have larger and more massive complexes as the 
scale for gravitational collapse (analogous to the Jeans 
length, but for turbulent rather than thermal velocities) 
is correspondingly larger. Complexes in HCG 31 were 
identified from the -B435 and composite images, visually 
and with the use of contour plotting, and photometered 
with rectangular apertures from the custom multidriz- 
zled images (1 pix = 0'.'03). For consistency with previ- 
ous work in the literature (see §3.6 and §3.7), no back- 
ground subtraction was done. The complexes were typ- 
ically larger than ~20 pixels in diameter (correspond- 
ing to ~150 pc), with outer boundaries ~10cr above the 
background. For the sake of comparison with other sys- 
tems, the focus was on identifying the largest and bright- 
est complexes (which are the most apparent both locally 
and at high redshift), and the sample is therefore incom- 
plete. There were 40 such complexes: 11 each in regions 
A+C, B, and G, 3 in region E, and 4 in region F, as 
shown in Figure 4 and 5 and listed in Table 2. The com- 
plexes in Region G correspond to the "optical knots" 
noted in ground-based unsharp-masked -B-band images 



by Verdes-Montenegro et al. (2005). A visual examina- 
tion of ground-based g and i images (Mendes de Oliveira 
et al. 2006) found no complexes of comparable size or 
brightness in the northern galaxies Q and R outside of 
the ACS field of view. 

The complex diameters in HCG 31 range from ^150 
to 875 pc. In all cases, they are significantly resolved, 
and while they may encompass several SC candidates, 
the SC candidate and complex sizes do not overlap. In 
nearby spirals, the sizes of the largest complexes scale 
with the absolute magnitude of the galaxy (Elmcgrecn 
et al. 1994, 1996) and range from ~T50 pc for Ms=-16 
to ^1500 pc for Mb = -21. The absolute V magnitude 
for HCG 31 A+C is -20.7 (for DM = 33.83; Conti 1991), 
with Mb 436 = -20.0 based on our HST image; its complex 
diameters range from 190-540 pc. B and G are measured 
to have Mb 435 of -18.5 and -18.9, with complexes ranging 
from 165-190 pc and 140-400 pc, respectively. These are 
comparable to the sizes of complexes in local spirals with 
similar absolute magnitudes. The complexes in the tidal 
debris of HCG 31 are anomalously large; region E, with 
Mb 435 = -16.5, has a complex with a diameter of 875 pc, 
and region F, with Mb 435 = -15.8, has a complex of 550 
pc. These values are comparable to what is observed in 
local Blue Compact Dwarf (BCD) galaxies and factors 
of 2-3 x larger than those in Sm and irregular galaxies of 
similar magnitudes (Elmegreen et al. 2006). 

2.4. Extended Sources 

To search for extended sources that might be iden- 
tified with either tidal dwarf galaxy candidates or dis- 
solved clusters, we ran the SExtractor software package 
(Bcrtin & Arnouts 1996) four times, each time with a 
different background smoothing size, which typically fa- 
vors objects of different sizes. After merging the catalogs, 
multiple detections were combined. The 'best' detection 
was chosen to be (in order) the roundest, biggest, and 
least stellar source. By visual inspection of the extended 
sources (particularly in the tidal debris of E, F, and H) 
identified by these means, no obvious dwarf galaxy can- 
didates were evident. The majority of the most signif- 
icant candidates can be identified cither with the com- 
plexes or marginally resolved SC candidates. Further- 
more, the -B435-V6O6 and ^06-^814 color distributions of 
the SExtractor sources were consistent with these pop- 
ulations. Therefore, we do not find a distinct population 
of extended sources to identify as either dwarf galaxy 
candidates or dissolved clusters. 

2.5. Evolutionary Track Models 

We used the evolutionary tracks of the Bruzual & 
Chariot (2003) instantaneous burst model assuming a 
Salpeter initial mass function (IMF) , a stellar mass range 
from 0.1-100 M Q , and 1/5 solar metallicity, consistent 
with the range of Z/Z Q =0.2-0.4 found by Lopez-Sanchez 
ct al. (2004) for this system. The spectra at each grid 
point in time have been convolved with our filter set to 
obtain colors and magnitudes. 

The Bruzual & Chariot (2003) models do not include a 
contribution from nebular emission. Emission lines such 
as Ha and [O III] can significantly affect the colors of star 
clusters with ages < 10 Myr. Given that the F606W fil- 
ter is quite broad, strong Ha emission at the redshift 
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of HCG 31 will boost the V<306 fluxes, thereby redden- 
ing the -B435-V6O6 colors and making the V606 _ ^8i4 colors 
bluer (e.g., Vacca & Conti 1992; Conti ct al. 1996). To 
calculate evolutionary tracks that incorporate a model 
of nebular emission, we used the Starburst99 19 (Lci- 
therer et al. 1999) emission-line strengths for Ha and H/3 
for the same IMF and metallicity burst as given above. 
Another strong line is [O III], which is not included in 
the Starburst99 models. To estimate the maximum 
expected contribution of [O III] to the colors, we used 
the highest value of 0.7 for the ratio of log([0 III] /H/3) 
observed in the KISS sample of low-mass star-forming 
galaxies (e.g., Salzer et al. 2005); these values have been 
included in the calculation of the nebular tracks. 



2.6. Intrinsic Reddening 

The colors of young star clusters are often reddened 
by the dusty clouds in which they are born. The mag- 
nitude and direction of an Aqqq—1 Galactic reddening 
vector (Cardclli ct al. 1989) is plotted in the color-color 
and color- magnitude diagrams (see Figs. 3 and 6). This 
may not be appropriate for these galaxies, given their 
lower mctallicities (Z/Z@ = 0.2-0.4; Lopez-Sanchez et al. 
2004) , but is conservative and consistent with convention 
in the literature. The V606 filter covers strong nebular 
emission lines such as [O III] and Ha from the youngest 
star clusters at the redshift of HCG 31, and as we showed 
in §2.5 the inclusion of this flux stretches the predicted 
evolutionary tracks down and to the left in Figure 3, 
nicely localizing very young star clusters still surrounded 
by interstellar material in color-color space. The redden- 
ing vector is nearly perpendicular to the tracks including 
nebular emission. Therefore, as long as extinction does 
not push too many star clusters below our selection limit 
of Mv 606 < —9, the youngest clusters will still be iden- 
tified as such. 20 This is not the case however, for older 
clusters which no longer ionize nearby interstellar gas. In 
this situation, shown as the solid line in the same Fig- 
ure, the reddening vector is nearly parallel to the tracks 
themselves, making it difficult to quantitatively estimate 
the total extinction on a cluster by cluster basis; to do 
so would require U -band data of comparable resolution 
(e.g., Whitmore et al. 1999). The colors therefore give 
upper limits for age estimates. 

We can get an independent assessment of extinction 
on larger physical scales from mid-infrared observations 
taken with the IRAC instrument on Spitzer. From the 
Spitzer IRAC images presented in Johnson et al. (2007), 
the strongest dust emission is present at the A+C in- 
terface. This is illustrated as the reddish-orange glow in 
Figure 1 (also see Fig. 4 of Johnson et al. 2007). At 8/zm, 
Galaxy B is quite faint, and the tidal debris is largely 
invisible; therefore we do not expect significant dust ex- 
tinction on average in those regions. The asymmetry 
in the disk of G is also evident in this representation. 
Therefore, in general we expect dust extinction will be 
most problematic in the main interaction region between 
A+C and the western star-forming ring of G. 

19 http:/ /www. stsci.edu/science/starburst99/ 

20 An assumption underlying this assertion is that ionizing stel- 
lar flux must be able to escape its natal dust cloud and still en- 
counter sufficient gas to create detectable nebular emission. 



2.7. Swift and GALEX UV Imaging 

We used UV observations taken with Swift /UVOT and 
GALEXto probe star formation activity in HCG 31 (see 
Tzanavaris et al., in preparation for details). Briefly, we 
calculated background-subtracted fluxes at the effective 
wavelengths of the U, UVW1, UVM2, UVW2 UVOT 
filters (3501, 2634, 2231, 2030A) and the NUV and FUV 
GALEX filters (2271 and 1528A) in the galaxies and 
tidal structures using the Johnson et al. (2007) apertures. 
These values have been corrected for Galactic extinction 
(Schlegel et al. 1998; Cardelli ct al. 1989) and are listed 
in Table 3. The errors include a systematic error that 
reflects the uncertainty in the zeropoints. No corrections 
for intrinsic extinction are included. 

Qualitatively, we note that the significant flux levels in 
the Swift UVW2 and GALEX FUV filters are indicative 
of stellar populations characterized by massive stars with 
ages of a few Myr. This is consistent with previous UV 
work that only covered a small area in the central region 
of A+C (Conti et al. 1996). The GALEX NUV data has 
been included as purple in Figure 1, and is consistent 
with the Swift UVM2 emission. We combine the UV 
and Spitzer 24/im fluxes to obtain star-formation rate 
estimates in §3.3. 

2.8. Multicolor Images 

Given the wealth of information from the optical, in- 
frared and ultraviolet images, we have constructed two 
full- field, multicolor images. The first, Figure 1, incor- 
porates images from five bands: GALEX (NUV), HST 
(-B435; V(306; aud ^814) aud Spitzer (8 /xm), and the sec- 
ond, Figure 2, is a three-color composite of HST data 
only. All data were converted to flux density units of erg 
s _1 cm~ 2 A -1 . These values were initially distributed 
to the greyscale range logarithmically. The Karma visu- 
alization package (Gooch 1996) was used to rotate the 
data and regrid the non-HST data to the HST pixel 
scale. With this, we also further adjusted the grey scales, 
for example, to emphasize the brightest regions in the 
non-HST data. The remaining data manipulation was 
done using GIMP. 21 Following the methods in Rector et al. 
(2007), color was assigned (-B435 = blue, V606 = yel- 
lowish green, I 8 i4 — red, NUV = purple, and 8 ^m = 
orange) and the images combined using the screen algo- 
rithm. Masks were applied to remove the ghost reflection 
and to reduce the UV emission from the foreground star. 
Selected regions of the; tail were brightened so that they 
would be visible, because requiring a dark background 
reduces the contrast of faint astronomical objects. 

3. RESULTS AND DISCUSSION 

3.1. SC Candidate Luminosity Functions 

The luminosity function of star clusters can be approx- 
imated by a power law, dN/dL oc L a , with a ~ —2, in 
many spirals, starbursts, and merging galaxies (see com- 
pilation in Whitmore 2003 and references therein). This 
strong similarity in the cluster populations across galax- 
ies types and with very different star formation rates led 
Whitmore (2003) to suggest that the basic properties, 

21 Gnu Image Manipulation Package (GIMP) written by Peter 
Mattis and Spencer Kimball and released under the GNU General 
Public License. 
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such as age, mass, and luminosity, of star cluster systems 
in nearby, actively star-forming galaxies could be broadly 
understood from simple statistics. In general, the lumi- 
nosity of the brightest young cluster (Larsen 2002; Whit- 
more 2003) in a galaxy scales with the total number of 
clusters in the system brighter than My= —9, regardless 
of whether the galaxy is quiescent or starbursting, dwarf 
or giant, interacting or isolated. Given a universal power- 
law luminosity function for young star clusters, galaxies 
with high star formation rates, and therefore larger star 
cluster populations, arc more likely to populate the lu- 
minous end of the luminosity function. 

In Figure 7, we construct the Vqog luminosity function 
for individual galaxies in HCG 31, a dynamical environ- 
ment with on-going, multiple interactions that is differ- 
ent than those studied previously. The power-law index 
a is calculated from the best-fit index /3 in these distri- 
butions as: a = 2.5 x /J + 1, and lies between « —1.8 
and —2.3 for each galaxy. To improve statistics, we also 
summed clusters from all galaxies together to obtain an 
index a « —2.1. This is well within the range found for 
cluster systems in very different galaxies. We show the 
relationship between the brightest cluster and the total 
number of clusters with My 606 < —9 for each galaxy in 
Figure 8. The line in this figure shows that the best 
fit relation determined by Whitmore ct al. (2007) from 
40 galaxies (not including any in compact groups) nicely 
describes the location of the HCG 31 galaxies in this pa- 
rameter space. The scatter in HCG 31 is also well within 
that observed. 

The discussion above assumes that dust has not signif- 
icantly affected the absolute magnitudes of the clusters, 
and so we consider its potential effects. The brightest 
clusters in each region are indicated with open stars in 
the color-color and color-magnitude diagrams for each 
region (Figs. 9 and 10). Their positions in color space 
(Fig. 9) are consistent with young (r < 100 Myr) or in- 
termediate (r=100-500 Myr) ages. If we assume that 
they are actually all very young (t < 10 Myr) and are in 
their current positions as a result of reddening, the in- 
ferred extinction would be A 6 o6^ 1 for all but E, whose 
colors can accommodate 1.5). Even in this worst 

case scenario, the absolute magnitudes of the most lumi- 
nous clusters would still be consistent with the Whitmore 
et al. (2007) relation, with the exception of the A+C clus- 
ter which would be about a magnitude brighter than the 
scatter. 

In sum, these results suggest that despite the unusual 
environment of a compact group, there are no obvious 
differences in the distribution of star cluster luminosities 
and numbers when compared with more typical environ- 
ments - which include spirals, starbursts, and mergers - 
locally. 

3.2. Ages and Spatial Distributions of Star Cluster 
Candidates 

The age distribution of star clusters provides impor- 
tant information on their formation and disruption, and 
can be determined by comparing the integrated cluster 
colors with the predictions of stellar population models. 
Star cluster candidates have been sorted according to the 
structures within HCG 31 they are spatially associated 
with as drawn in Figure 2. From inspection, the ma- 
jority of the SC candidates are spatially coincident with 



galaxies and tidal features in HCG 31, with only 4.4% 
(19 of 424) outside of these boundaries. Furthermore, 
the SC candidates positions within the color-color plot 
(Figs. 3 and 9) are generally consistent with instanta- 
neous burst evolutionary tracks. Notably, a large fraction 
of the brightest clusters in each region are in the lower left 
section of Figure 3 as expected for the youngest clusters 
(r < 10 7 yr) and typical of galaxies currently forming 
stars. The few cluster candidates which fall up and to 
the left of the region where the nebular emission tracks 
connect to the rest of the models in the Figure 3 are 
likely to be < 10 7 yr, blue clusters with little nebular 
emission - clearly a minority of the youngest clusters. 
Though galaxy A+C is clearly where the largest frac- 
tion of young clusters are found, galaxies B and G also 
have significant populations of young clusters with colors 
consistent with nebular emission (see Fig. 9). 

Many of the cluster candidates in the HCG 31 galaxies 
fall in a clump in color-color space around V(306~^8i4~ 0.4 
and £?435-V(506~ 0-2, near the region where the tracks 
which include nebular emission join those which do not. 
Given the typical photometric errors, the colors of these 
objects are consistent with ages of 10-500 Myr. In addi- 
tion, young clusters often have at least modest amounts 
of extinction associated with them or their parent galax- 
ies, and so younger clusters (without nebular emission) 
can be pushed into this region. Unfortunately, we are 
unable to correct for extinction on a cluster by cluster 
basis because dust reddening is degenerate with age in 
the filter combination used here. Therefore, we can place 
only an upper limit on the estimated age of any given 
cluster that falls along the non-nebular track. Despite 
this limitation, the color-color diagrams reveal interest- 
ing features. 

Our cluster sample consists of all objects with 
Mv 606 < — 9, where we have only corrected for fore- 
ground extinction. Whitmore et al. (1999) have shown 
that selecting only objects with My < —9 returns a very 
clean sample of star clusters, because this limit is brighter 
than all but the most luminous stars. Because clusters 
fade rapidly with age, a number of them will naturally 
fade below a given magnitude or luminosity and be lost 
from the sample. Therefore, this sample is not sensitive 
to lower mass clusters at older ages (though see §3.4 for 
a discussion of globular cluster candidates), which can 
partially explain the lack of many clusters with redder 
colors. 

We can interpret the color (or equivalently, a crude 
age) distribution of clusters in Figure 3 by consider- 
ing the results from simple statistical cluster popula- 
tion models (e.g., Whitmore et al. 2007). Specifically, 
the formation and evolution of star cluster systems can 
be approximated statistically by following simple formu- 
lae, because star cluster systems in galaxies of different 
types show regularities, such as very similar power-law 
luminosity functions as discussed above. Chandar, Whit- 
more, & Fall (in preparation) predict the number of clus- 
ters expected for a simple model which assumes that clus- 
ters form at a constant rate with an initial power-law 
mass function, dN/dM oc M~ 2 . Following this simple 
prescription, if no clusters are destroyed and assuming a 
single V(so6 luminosity limit to mimic our selection proce- 
dure, the number of clusters increases each decade in age, 
by factors of w 2 (10 7 - 10 8 yr) and w 4 (10 8 - 10 9 yr) 
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over the number of 10 6 — 10 7 yr clusters. Even taking into 
account the dimming of clusters with time, the number 
of detectable older clusters in bins of 1 dex in log(r) is 
still expected to grow if no clusters are destroyed. The 
distribution of cluster colors in HCG 31 galaxies how- 
ever, shows that the opposite trend is observed, i.e., the 
number of clusters falls off rapidly with color (age) for 
HCG 31 galaxies redward of the clump of clusters men- 
tioned above at estimated ages wfewxlO 7 yr. Therefore, 
a scenario with no cluster destruction is not viable unless 
the star formation rate has increased significantly in the 
very recent past. 

This rapid decline in the number of clusters with red- 
der colors, observed for all of the galaxies, results from 
a changing history of formation or ongoing disruption of 
the star clusters, or to some combination of both. A rapid 
decline in the number of clusters observed with age (or 
color) has been observed previously in several galaxies, 
including the Antennae (e.g., Fall et al. 2005; Whitmore 
et al. 2005), the solar neighborhood (Lada & Lada 2003), 
and the Small Magellanic Cloud (Chandar et al. 2006). 
The early, rapid disruption of clusters at ages r < 10 8 yr 
is plausibly dominated by two physical processes internal 
to the clusters themselves: the disruption of star clusters 
by the expulsion of interstellar material left over from 
star formation, and mass-loss from continued stellar evo- 
lution (Fall, Chandar, & Whitmore, submitted). Because 
these internal processes are largely independent of envi- 
ronment, they should disrupt clusters in a similar way 
in most galaxies. This picture predicts that the number 
of clusters in a luminosity-limited sample should decline 
with age (or towards redder colors), very similar to the 
pattern observed here for several galaxies in HCG 31, 
even with a constant star formation rate. Therefore, de- 
spite the seemingly unusual environment of a compact 
group, there are no obvious differences in the distribution 
of star cluster colors, luminosities, and numbers when 
compared with more typical environments. 

An alternative explanation is that all of these galaxies 
in HCG 31 have experienced a very large, recent, coordi- 
nated burst of cluster formation. These galaxies have a 
range of star formation rates, (see §3.3), that appear to 
reflect their interaction history. For example, the most 
actively star-forming galaxies, A and C, are experiencing 
a direct interaction. Though this interacting/merging 
galaxy pair is forming stars at a higher rate than galaxy 
G, and these galaxies do not resemble each other mor- 
phologically at all, the general distribution of the clusters 
colors (ages) is quite similar. In a typical environment, 
there is no reasonable physical mechanism for coordinat- 
ing star formation on such large scales given a character- 
istic speed of bulk motions in the ISM of ~ 100 km s -1 . 
However, in HCG 31, A+C, G, and B as well as the 
tidal features E, F, and H, are connected by a common 
H I envelope of tails and bridges, and each region is un- 
dergoing a significant rate of star formation given their 
galaxy types (see §3.3). Therefore, while the SC can- 
didate distributions of colors and luminosities are likely 
to be significantly shaped by the mechanisms of creation 
and disruption common to all star-forming galaxies, it is 
clearly the case that star formation in this group as a 
whole has been recently elevated as a result of the mul- 
tiple interactions. 



3.3. Star Formation Rates 

Massive, young stars emit light at UV wavelengths thus 
offering a direct measure of ongoing star formation. Some 
fraction of this light is readily absorbed by dust and then 
re-emitted in the IR. Therefore, UV emission lost to the 
observer can be recovered using IR data. Following the 
methodology of Tzanavaris et al. (in preparation), we 
estimate a total star-formation rate (SFR) as the sum 
of unobscured (UV-derived) and obscured (IR-derived) 
components. To the Spitzer 24 ^.m photometry for the 
HCG 31 galaxies published in Johnson et al. (2007), we 
applied the calibration of Rieke et al. (2009) to obtain the 
fraction of star-formation rate (SFR) probed by the IR. 
Next, we applied the Kennicutt (1998) calibration to the 
Swift /UVOT photometry to obtain the dust- unobscured 
fraction of the SFR. The sum of these two gives a total 
SFR effectively corrected for intrinsic dust extinction. 
Our estimates are listed in Table 3. 

The total UV+IR-derivcd SFR for all galaxies is ~ 
10.6 M Q yr -1 , with the largest contribution (> 75%) 
due to the combined emission from A+C and E (unre- 
solved in the 24 /zm Spitzer images). This is in reason- 
able agreement with previous estimates using different 
SFR indicators, such as Ha, far-IR, 60^m and 1.4 GHz, 
which give estimates ranging from 2-8 Mq yr -1 (Dopita 
et al. 2002; Richer et al. 2003; Lopez-Sanchez et al. 2004; 
Lopez-Sanchez & Esteban 2007) for the whole group. 

As an independent means of constraining the SFR with 
the new HST data, we examined the SC candidates with 
colors consistent with < 5 Myr star clusters with nebular 
emission (see §2.1 and the dashed green line in Fig. 3). 
Summing the masses of these clusters and dividing by 
the integrated time on the nebular evolutionary track 
will give an estimate of the average cluster SFR over this 
epoch. To convert the colors to an age, the -B435-V6O6 vs. 
V606~^8i4 axes were rotated so that the reddening vector 
was vertical. For most of the track, the horizontal posi- 
tion in this rotated color space increases monotonically 
with time. The colors of the SC candidates designated 
as nebular were then projected onto the nebular evolu- 
tionary track, and the distance to the track was recorded 
as the intrinsic reddening of the cluster. (Negative dis- 
tances were set to 0.) This gave a mean extinction value 
for the sample of v4 6 06 = 0.74. The Ve;o6 magnitude of 
each SC candidate was then corrected for its intrinsic 
reddening and converted to a mass using the models (see 
the dashed curves in Fig. 6). Following this algorithm, 
we obtained a cluster SFR of <~ 1 M yr -1 , approxi- 
mately an order of magnitude lower than the estimate 
from the completely independent UV+IR emission. If 
interpreted literally, this suggests that, given our incom- 
pleteness to low mass and highly cxtinctcd clusters, at 
least one-tenth of newly formed stars are still associated 
with their natal star clusters at these early times. A 
caveat to this statement is the unknown systematic un- 
certainty in translating color and magnitude to age and 
mass. 

3.4. Globular Clusters in HCG 31 

Some of the point sources from Figure 3 are expected 
to be old (r > 10 Gyr) globular clusters (GCs) in the 
HCG 31 system. Though GCs have been studied in only 
a few HCGs to date (Barkhouse et al. 2001; Da Rocha 
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et al. 2002b) , they can be used to provide information on 
the formation and early evolution of such systems (e.g., 
Brodie & Strader 2006). 

Use of color-color diagrams to identify such clusters 
has proven very effective (e.g., Rhode & Zepf 2001, 
2004) . The color-color diagram for all point sources with 
V606< 26 (or My 606 < —7.8 where we expect over 50% 
completeness for old clusters) in the entire HCG 31 field 
is plotted in Figure 11. We have relaxed the magni- 
tude limit to Mv 606 < —7.8 because for red objects we 
are not concerned that single stars could be this bright, 
in contrast to the discussion of the young, blue clusters 
in §3.1 and §3.2. To define the color range for globular 
clusters, we have used the reddening corrected B — V 
and V — I colors for 97 Milky Way globular clusters 
(those with BVI magnitudes) from the Harris (1996) 
catalog. The color range 0.55 < (B — V)q < 1.05 and a 
0.2 mag wide swath in (V — I)o centered on the relation 
(V - I) = 0.82(B - V) + 0.36 was found to contain 
94% of the clusters. This GC bounding region was then 
transformed to our ACS filter set (Sirianni et al. 2005) 
and is shown in Figure 11. 

A total of 43 point sources with Veo6< 26. (or My 606 < 
—7.8 at the distance of HCG 31) have error ellipses that 
overlap the bounding region, and are thus considered 
candidate globular clusters in HCG 31. The vast ma- 
jority of these objects are centered on the primary galax- 
ies in the group. In A+C, B, and G, they are spread 
throughout the galaxies, and so the majority of them are 
not expected to be simply younger, reddened clusters. 
The latter would be concentrated at the dust emission 
peaks at the A+C interface and along the western arm 
of G (see §2.6 and Fig. 1). Although some background 
elliptical galaxies and foreground Milky Way dwarf stars 
with similar colors are expected to contaminate our sam- 
ple, the paucity of such objects over most of the field 
(coded as 'Other' in Figure 11) shows that the vast ma- 
jority of objects are likely true GCs in HCG 31. We have 
used the Besangon Milky Way stellar population model 
(Robin et al. 2003) to estimate the number of foreground 
Milky Way stars with magnitudes and colors consistent 
with our GC selection region; we expect a total of 2 to 
3 objects over the entire ACS field for the stellar density 
of 0.23 objects arcmin -2 , consistent with the number of 
objects (2) that already clearly lie outside the defined 
HCG 31 regions. Therefore, background contamination 
is not significant in our candidate GC sample. Inspection 
of the color-color plot of strictly defined point sources 
supports our conclusion that the vast majority of ob- 
jects with nebular colors are marginally resolved; only 
two point sources with nebular colors (as defined in Fig- 
ure 3) and My 606 < —9 satisfy the point source selection 
criteria. 

Studies of GCs in late-type galaxies such as these are 
made extremely difficult by the rapidly varying back- 
ground close to the centers of the galaxies, and there is a 
real possibility that some objects suffer additional inter- 
nal extinction. As a result, some of our GC candidates 
may well be younger objects that have been significantly 
reddened (note the reddening vector in Figure 11). Fur- 
thermore, the detection of objects close to the galaxy 
centers may be adversely affected by photometric incom- 
pleteness - the 90% completeness level for point sources 
within the HCG 31 subsections (see §2.2) are significantly 



brighter than outside the galaxies. 

From inspection of the color-magnitude diagrams for 
the GC candidates, we find that while many of the can- 
didates are spatially close to their parent galaxies, the 
large number close to the average 90% completeness lim- 
its rather than the brighter galaxy core limits, indicates 
that the former are more appropriate. As a result, we 
have adopted a fraction / = 0.9 ±0.1 for all objects with 
V(506< 26. This is certainly not applicable for a small 
number of objects very close to the galaxy centers, but 
these will not affect the basic results from our study. 

A further test on the effects of incompleteness can be 
gleaned from the Vqog luminosity function of the candi- 
dates, shown in Figure 12. While the number of GC can- 
didates is small, the luminosity function does show the 
general property expected for a globular cluster system: 
a monotonic rise towards the turnover of the GC luminos- 
ity function. Also plotted is the expected GC luminosity 
function assuming a peak at My ~ —7.4 (Vqoq = 26.3 
for DM = 33.83, V - I = 1.0 ± 0.3), and a dispersion 
a = 1.2 (Harris 2001). That the data continue to rise 
very sharply to Veoe= 26 again suggests that photomet- 
ric incompleteness due to varying galactic backgrounds 
are not strongly affecting our numbers. 

From the expected GC luminosity function (as plot- 
ted in Figure 12), we estimate that the GCs brighter 
than ^606= 26 are 40 ± 7% of the total GC distribu- 
tion, where the errors reflect added uncertainties of 0.2 
mag in the GC luminosity function dispersion, and an 
uncertainty of 0.2 mag in the location of the peak mag- 
nitude. To estimate the total number of GCs in each 
galaxy in HCG 31, we first correct our observed num- 
bers N bs (and associated Poisson errors) by our bulk 
completeness fraction / = 0.9 ±0.1, and then correct for 
the fraction of the luminosity function we can observe. 
The final, corrected total number of GC candidates are 
given in Table 1, where we have also included the single 
GC candidates superposed on the tidal features E and 
F. We have made no corrections for the (small) back- 
ground contamination derived above. Fortunately, both 
the large distance to HCG 31 and the compactness of the 
GC systems in these low-mass galaxies mean we are es- 
sentially observing the entire GC system, and so we make 
no corrections for any clusters outside the ACS field of 
view. 

Column 6 of Table 1 shows the specific frequency Sn, 
or the number of GCs normalized to the U-band lumi- 
nosity M v = -15 (Harris & van den Bergh 1981). 22 The 
values of Sn for each of the larger galaxies A+C, B and 
G are all consistent with the lower values expected for 
late- type galaxies in other environments (Sn ~ 0.5 — 0.9; 
e.g. Goudfrooij et al. 2003; Chandar et al. 2004; Rhode 
et al. 2007). From this, there is no evidence that these 
galaxies have undergone previous mergers that formed a 
significant population of new massive star cluster popu- 
lations. 

3.5. Intragroup Stellar Populations? 

While the number of GCs per unit mass (the T parameter 
from Zepf & Ashman (1993)) is now commonly used to compare 
numbers of clusters between differing galaxy types, all galaxies in 
HCG 31 are late-type galaxies and the 5jv values can be compared 
directly. 
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The lack of a significant number of GC candidates 
far outside the galaxies may also indicate a dynamically 
young system, because galaxy interactions are expected 
to liberate some stars (and star clusters) from the gravi- 
tational potentials of their parent galaxies. Indeed, some 
compact groups do show significant amounts of diffuse 
light (e.g., Pildis et al. 1995; White et al. 2003; Da Rocha 
& Mendes de Oliveira 2005; Da Rocha et al. 2008). Does 
a significant diffuse component (in the form of intergalac- 
tic GCs) exist in HCG 31? As there are only two objects 
that lie significantly outside the small globular cluster 
systems concentrated in the primary galaxies (broadly 
consistent with the expected number of foreground Milky 
Way stars), there is no significant population of intra- 
group GCs present in our field. However, we can put a 
meaningful upper limit on the presence of such a popu- 
lation. 

To quantify this further, we first assume that both GC 
candidates outside the galaxies are indeed intragroup sys- 
tems. We define the fraction of intcrgalactic cluster light 
as fw = L IC /{L ga i + L IC ), where L gal is the combined 
y-band luminosity of the HCG 31 galaxies and Lie is 
the total luminosity of the intergalactic GC population. 
Using the My value of Conti (1991) converted to our 
value of the DM , we find a combined My g roU p = -20.7, 
or Ly ~ 15.5 x 10 9 L Qi y. Given the two intergalactic 
GC candidates located far from the HCG galaxies, and 
subtracting the expected contamination of two stellar ob- 
jects, we have an observed Niac,obs consistent with zero, 
or (correcting for the missing part of the luminosity func- 
tion) an upper limit of 6 intergalactic GC candidates in 
the ACS field. Adopting a mean group specific frequency 
Sm = 0.6 ± 0.3 (appropriate for the late type galaxies in 
HCG 31), this upper limit corresponds to a F-band lumi- 
nosity of ~ 0.9 x 10 9 i Qj y. Thus we place an upper limit 
of 5% on the intergalactic light in HCG 31. This low 
value of the intergalactic light fraction, when combined 
with previous work (where intragroup light fractions in 
other HCGs can be as high as almost 50%) is suggestive 
that the range of galactic properties, physical conditions, 
and dynamical history in the compact group environment 
may be extremely important in understanding the differ- 
ent processes that shape group evolution. 

While our estimate above suggests there is very little (if 
any) intragroup light in HCG 31, a caveat here is that all 
galaxies in the group are rather low mass, with compact 
GC systems concentrated over a small (r < 10 kpc) radial 
extent (Rhode et al. 2007), so while it is possible that the 
interaction history of HCG 31 is so recent as to not have 
liberated much of an intergalactic population (consistent 
with the other observations presented here), it is unclear 
whether many GCs would be readily removed from such 
spatially concentrated GC systems in any case. 

3.6. Large-Scale Star Formation 

Star-forming complexes represent the largest units of 
star formation in a galaxy (see review by Efrcmov 1995); 
they are typically kpc-scale regions encompassing sev- 
eral clusters. In local spirals, their masses and sizes 
approximately follow a power law distribution, and the 
largest ones scale with the host galaxy's absolute magni- 
tude (Elmegrecn et al. 1994, 1996). Complexes form from 
large-scale gravitational instabilities that depend on the 
local velocity dispersion, whether the region forms in a 



disk or a tidal stream. In disk galaxies, the Jeans length 
is about equal to the disk thickness, and complexes often 
appear as rather regularly spaced beads of star forma- 
tion in normal spirals. Because a galaxy's largest com- 
plexes can be resolved even in high redshift galaxies, they 
are useful probes for comparisons of large-scale star for- 
mation in a variety of environments. In order to elimi- 
nate possible size-of-sample effects (discussed by Selman 
& Melnick 2000 but disputed by Vicari et al. 2002 and 
Elmegreen et al. 1994), complexes should be compared 
among galaxies with comparable absolute magnitudes. 

Figure 10 shows a color-magnitude diagram includ- 
ing the HCG 31 complexes, coded according to their re- 
gions. Superposed are the evolutionary tracks discussed 
in Section §2.5. The complexes in region F are too blue 
in V606~^8i4 f° r the basic tracks, and require the mod- 
els that include nebular emission. We have used these 
tracks for all of the complexes, because a comparison 
with the Ha image of HCG 31 (Johnson & Conti 2000) 
shows that all of the complexes are associated with emis- 
sion. Many of the complexes lie within the Ha knots 
identified from ground-based images by Iglesias-Paramo 
& Vilchez (1997). From their positions in the color- 
magnitude diagram, the inferred complex masses range 
from ~ 5 x 10 5 to 2 x 10 7 M Q ; overall, the complexes aver- 
age about 1O 6 M0. The most massive ones are in E and in 
A+C; two of them in the central regions of A+C prob- 
ably encompass more than one complex blended along 
the line-of-sight. The complexes are young, as expected 
by their Ha emission, with log(ages) ranging from about 
6.5 to 7.5 years. The oldest ones are in regions A+C and 
B, in the central parts of those galaxies, which inference 
may be due partly to blending with an older background. 

The luminosities, L, of individual complexes scale ap- 
proximately with their diameters, D, as L oc D 18 . This 
is essentially the same relation as that found for com- 
plexes in nearby quiescent spirals (Elmegreen et al. 2006; 
Elmegreen & Salzer 1999) and also in the gently interact- 
ing galaxy pair IC2163/NGC2207 (Elmegreen et al. 2001) 
at a distance of 35 Mpc. However, the absolute magni- 
tudes of the complexes are ~2 mag brighter in HCG 31 
than in IC 2163/NGC 2207 (even though those galaxies 
are 1 and 2 magnitudes brighter than HCG 31 A+C, re- 
spectively) and in quiescent spirals with comparable or 
brighter total absolute magnitudes. The masses of the 
complexes in HCG 31 therefore are also higher, by a fac- 
tor of ^5-10. This suggests that the velocity dispersions 
may be higher in HCG 31 than in IC 2163/NGC 2207. 
The complex masses in HCG 31 resemble those in blue 
compact dwarf galaxies, which are also larger than those 
in normal spiral galaxies with the same absolute magni- 
tudes. This result was inferred to be a consequence of 
higher disk turbulence in blue compact dwarf galaxies 
(Elmegreen et al. 1996). 

Ha Fabry-Perot observations by Amram et al. (2007) 
show velocity dispersions of typically 15-20 kms -1 over 
most of the regions in HCG 31, with values up to 
30kms _1 in the vicinity of the complexes. Elevated 
values in the immediate vicinity of the complexes prob- 
ably result from extra turbulence generated by associ- 
ated energetic events such as stellar winds and super- 
novae. Maps by Richer et al. (2003) with lower spatial 
and spectral resolution show velocity dispersions of 45 to 
95 km s _1 over the main disk of galaxies A+C, with fairly 
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constant values of ~50kms _1 in regions B, E, and F, and 
slightly higher values in region G. Although Ha veloc- 
ity dispersions do not necessarily reflect the kinematics 
of the neutral gas out of which the complexes formed, 
the observed high large-scale Ha velocity dispersions in 
HCG 31 suggest that the disks and the tidal debris are 
more turbulent than the 5-10 km s^ 1 typical of quiescent 
spiral galaxy disks; this is consistent with the presence 
of more massive complexes. 

Star formation is commonly observed in tidal features 
of interacting galaxies (e.g., Schweizer 1978; Mirabel 
et al. 1992). Numerical simulations indicate that 
complexes in such systems can form from interaction- 
triggered gravitational instabilities in the gaseous compo- 
nent (e.g., Elmegreen et al. 1993; Wetzstein et al. 2007). 
Complexes observed in tidal tails in the local interacting 
systems NGC 4485/90 (Elmegreen et al. 1998), the Leo 
Triplet galaxy NGC 3628 (Chromey et al. 1998), and the 
Tadpole galaxy UGC10214 (Tran et al. 2003) all have 
masses of ~ 10 6 M Q , similar to those in HCG 31. How- 
ever, all of these galaxies are 0.6 to 1.3 mag brighter than 
A+C, the brightest of the HCG 31 galaxies, and so the 
HCG 31 complexes are more massive relative to their 
galaxy luminosity. 

Another exceptional interacting system is Arp 285, 
which consists of two galaxies with a tidal tail containing 
beads of star-forming complexes observed with Spitzer 
(Smith et al. 2008). These galaxies are each ~ 0.8 mag- 
nitudes fainter than A+C; their tail clump masses are 
< 10 6 M Q , but they also contain a bright spot in the 
disk with 10 8 M . 

3.7. Comparison with Complexes in High Redshift 
Galaxies 

Complexes have been observed with HST in high 
redshift galaxies in the GEMS and GOODS fields 
(Elmegreen & Elmegreen 2006; Elmegreen et al. 2007a) 
and in the Ultra Deep Field (UDF; Elmegreen & 
Elmegreen 2005; Elmegreen et al. 2005, 2007b, 2009), 
with sizes in the kpc range. In the standard ACDM cos- 
mology (Spergel et al. 2007), 1 kpc corresponds to ~ C/13 
at z = 1, and has a nearly constant angular size over the 
redshift range 0.7 < z < 4. At closer redshifts, the angu- 
lar size decreases approximately linearly with increasing 
redshift. With 0'.'03 pixels in the ACS images, a kpc- 
size region corresponds to 3-4 pixels, so complexes much 
smaller than this could not be measured in the UDF 
galaxies. On the other hand, if there were complexes sig- 
nificantly larger than this in nearby galaxies, they would 
be easily discernible. In other words, the largest com- 
plexes can be identified and compared in both local and 
high redshift galaxies. The high z galaxy complexes have 
masses typically 10 7 to a few 1O 8 M from z = 1 — 4, de- 
creasing to about 10 6 M Q at z = 0.1 — 0.2, as shown in 
detail in figures in the previous papers. For the high 
z cases, masses are estimated by comparing the com- 
plex position in color-magnitude space (cf. Fig. 3) and 
comparing with evolutionary tracks that span a range of 
reasonable star formation decay rates. In HCG 31, the 
masses are estimated by comparison with single popu- 
lation, instantaneous burst evolutionary models. Given 
the colors of the HCG 31 complexes, the inferred masses 
are not highly sensitive to this assumption which results 
in a factor of approximately two uncertainty (cf. Fig. 10, 



Elmegreen et al. 2007a). 

Among obviously interacting galaxies, which include 
M51-type galaxies as well as strongly interacting, gas- 
rich pairs with long tidal arms (Elmegreen et al. 2007a), 
the complex masses are 10 5 to 10 6 5 M Q at nearby red- 
shifts up to z = 0.125, and 10 7 - 1O 9 M at redshifts 
1 < z < 1.4. Faint tidal features are difficult to discern 
at higher redshifts because of cosmological dimming. 

In the GEMS and GOODS samples, galaxies at z < 
0.2 mostly have rest-frame My fainter than -18, because 
the small volume covered to that depth does not include 
many intrinsically brighter galaxies; beyond z = 0.2, the 
galaxies measured have rest-frame My from -18 to -22. 
As in local galaxies, the largest complexes are dependent 
on the absolute (rest-frame) magnitudes of the galaxies, 
so the main reason the lower redshift GEMS and GOODS 
complexes are less massive is because their host galaxies 
are intrinsically fainter. The largest complex masses have 
a spread of about a factor of 10 for a given absolute rest- 
frame magnitude, ranging from ~ 10 6 5 M Q for rest-frame 
My of -16 to ~ 10 8 M for -20.5 mag (Elmegreen et al., 
in prep.). These masses are a factor of ^10 larger for 
a given absolute magnitude than for complexes in local 
galaxies. In order to compare the complexes in HCG 31 
with those in high redshift galaxies, a Gaussian blur of 30 
pixels was applied to the HCG 31 images, corresponding 
to pixels covering about 250 pc. This scale simulates its 
appearance if at 1 < z < 4. 

The HCG 31 complexes were re- measured based on 
their boundaries in the blurred images. Some complexes 
were not recognizable in the blurred images because they 
were too faint or too small; a few were blended with other 
complexes that were too close. Most of the lowest mass 
complexes disappear, although the remaining complexes 
still range from ~ 10 5 to 10 7 3 M Q . A comparison of 
the full resolution and blurred complexes shows that the 
blurred ones typically are only ^20% more massive, even 
when they contain blended complexes; the colors differ 
by at most ~ 0.2 mag, so the inferred ages are about the 
same. The diameters range from about 300 to 650 pc 
(in one case, the diameter measured was smaller than in 
the unblurred image, due to uncertainty in placing the 
boundary.) These complexes arc similar in mass and size 
to those in the interacting GEMS and GOODS galaxies 
of comparable rest-frame magnitudes. 

The larger masses in the high redshift galaxies are pre- 
sumed to be the result of high turbulence in the disks; 
for example, kinematics inferred from spectroscopic ob- 
servations of a z — 1.6 galaxy, UDF6462, indicates ve- 
locity dispersions of 40kms _1 (Bournaud et al. 2008), 
while numerical simulations reproducing the large-scale 
star formation (Bournaud et al. 2007) also require values 
several times higher than in local quiescent galaxies. 

Thus, the complexes in HCG 31 can be understood in 
terms of large-scale star formation from gravitational in- 
stabilities in the gas driven by turbulence. Their sizes 
are comparable to what is expected based on complexes 
in nearby spirals. Their masses are similar to those in 
z ~ 1 galaxies of comparable magnitude, and are in- 
termediate between the smaller ones found in gently in- 
teracting nearby systems and the larger ones found in 
presumably more turbulent high redshift counterparts. 

3.8. Underlying Stellar Populations 
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An outstanding issue from previous HST studies is 
whether F contains an older stellar population in ad- 
dition to the current generation of young stars. From 
an examination of the region-specific color-color plot 
(Fig. 9), none of the star cluster candidates or complexes 
in F have colors consistent with > 1 Gyr populations, 
and there is only one globular cluster candidate within 
the larger F boundary. In general, it is hard to see in- 
dividual, old stars directly (with the exception of GCs), 
due to their faintness. The diffuse light within galaxies 
offers an opportunity to study fainter stellar populations 
than those described in the previous sections and can 
give us a handle on the population of older stars. 

To search for a distributed, underlying, older stellar 
population, we examined a processed -B435-/814 image of 
the entire ACS field of view. The image, presented in 
Figure 13, was created from the point-source subtracted 
images generated by the SC candidate PSF-fitting pho- 
tometry. It has been Vega-magnitude calibrated and 
boxcar-smoothed (with a 13 pix [180 pc] window) and 
corrected for Galactic extinction. We use -B435-/814 to 
remove the confusion of nebular emission prominent in 
V(;o6) and also because i?435-/si4^ 1-0 is a good indica- 
tor of stellar populations with ages > 500 Myr, barring 
intrinsic reddening (see Figure 13). 

In actively star forming galaxies, the diffuse light will 
be a mix of different stellar populations. We might ex- 
pect the light to be dominated by the youngest (and also 
brightest) stars that are present. Therefore, we will typ- 
ically be able to put a lower limit on the age of stars 
contributing to the diffuse light in the galaxies. 

From inspection, diffuse tidal debris with -B435- 
^814^0.4-0.8 apparently pulled out from A+C is evident 
to the north, as is some debris stretching to the south and 
overlapping with E and H. These colors are consistent 
with stellar ages > 100 Myr in this debris. The overall 
colors of A+C and B are in the range -8435-/814=0.4- 
1.0, with bluer (S435-7gi4=-0.1 to 0.4) regions coinci- 
dent with the locations of the youngest SC candidates. 
While the inner disk of G is blue (-6435-/814^ 0.4), the 
location of the complexes to the west stand out with 
i?435-/8i4~-0.1 to 0.4. The outer disk of G is signifi- 
cantly redder than the other structures in the group, with 
/?435 _ ^8i4=0.7-1.2. From the -B435-/814 color-magnitude 
diagram (Fig. 14) this indicates that the intermediate age 
range in A+C, B, and G is at least several 100 Myr, up 
to ~ 1 Gyr. The presence of GC candidates (§3.4) also 
supports the presence of ancient stars in these galaxies. 

Notably, F shows no evident diffuse emission consis- 
tent with an older stellar population; the reddest regions 
have i?435-/8i4~ 0.2-0.3, consistent with stars < 10 Myr. 
Given the low surface brightness of this material, it can- 
not be hiding a substantial population of old stars. While 
we cannot rule out that there is some low level stellar 
population with ages of 100 Myr and older, these ob- 
servations show no evidence for such a population, and 
probably give the best evidence to date that F started 
forming stars only recently. This implies that the ma- 
terial comprising F is almost entirely 'fresh', and F is 
overall a very young structure. 

4. SUMMARY AND CONCLUSIONS 
4.1. Modes of Star Formation 



An in-depth examination of the colors of structures at 
all scales of star formation, from compact star cluster 
candidates, to large star-forming clusters, to the diffuse 
underlying emission, indicates that the entire HCG 31 
system is suffused with recent and ongoing star forma- 
tion. Star cluster candidates with nebular colors arc 
present throughout HCG 31, specifically concentrated 
in the inner disk and western ring of G, through the 
midplane of B, in the interaction region of A+C, and 
across E, H, and F. While the prominence of SC can- 
didates with nebular colors is not sufficient evidence of 
a recent enhancement of star formation activity, other 
available information points in that direction. In par- 
ticular, while the SFR is enhanced relative to galaxies 
with similar masses, the sum total of H 1 in the sys- 
tem is not noticeably depleted. Crudely, this sets an 
upper limit to the timescale of enhanced star formation 
of ~ 300 Myr, as in this amount of time a significant 
fraction (~ 20%) of the H 1 would have been converted 
to stars. Furthermore, the tight constraints on intra- 
group GC candidates and diffuse light, and the consis- 
tency of the GC luminosity distribution with the low end 
of that expected for similar galaxies argue against pre- 
vious epochs of interaction-triggered star formation. In 
simulations of mergers of pairs of bulgeless disk galaxies, 
the typical timescale for the SFR to be elevated by an 
order of magnitude or so is ~ 150 Myr, and the most in- 
tense burst is triggered shortly after the first encounter 
(Mihos & Hcrnquist 1996). If these results are applica- 
ble to the more complicated dynamics of the compact 
group environment, then the ongoing burst may be of 
quite recent vintage. The current high rate of star for- 
mation is not sustainable long term (at most a Gyr or 
so) given the finite reservoir of cold gas of 1.6 x 10 10 M Q 
(Verdes-Montenegro et al. 2005). 

On the scales of star clusters, both young and globu- 
lar, the HCG 31 system is consistent with local galaxies 
of similar star formation rates and masses, as shown by 
the luminosity functions (Figs. 7 and 12). Furthermore, 
the values of brightest cluster My eoe as a function of 
number for distinct regions in HCG 31 match those of 
galaxies in the local Universe with a range of properties 
(Fig. 8). This supports the claim that star cluster for- 
mation (and destruction) follow universal patterns that 
scale with star formation rate, but are otherwise insensi- 
tive to their large-scale environment. 

However, the SC candidate populations of all com- 
pact groups studied to date are not the same. Stephan's 
Quintet (HCG 92) is another Hickson Compact Group 
with active star formation as a result of strong gravita- 
tional interactions (e.g., Xu et al. 1999; Sulentic et al. 
2001). Specifically, the galaxy NGC 7318b is blueshifted 
by 900 kms -1 with respect to the rest of the group galax- 
ies, and is currently ploughing through the neutral intra- 
group medium (Shostak et al. 1984). This has triggered a 
large-scale (~ 40 kpc) ridge of star formation coincident 
with a radio and X-ray bright shock (e.g., van der Hulst 
& Rots 1981; Pietsch et al. 1997; Trinchieri et al. 2003). 
With HST WFPC2 three-color imaging, Gallagher et al. 
(2001) identified a population of young, luminous SC 
candidates that traces the large-scale shock. In contrast 
to HCG 31, however, none of these SC candidates has 
colors consistent with strong nebular emission, though 
clusters apparently as young as ~ 2-3 Myr are seen. Af- 
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ter accounting for the brighter absolute magnitude limit 
and different filters of the WFPC2 observations, SC can- 
didates in Stephan's Quintet with comparable colors and 
absolute magnitudes (of My < —10) to those in HCG 31 
would have been detected if present. There are in total 
~ 30 star cluster candidates in Stephan's Quintet with 
comparably young ages (t < 10 Myr), the vast majority 
of which lie along the large-scale shock. From Spitzer 
spectroscopy of the shock ridge, strong, broad H 2 emis- 
sion lines, weak polycyclic aromatic hydrocarbon emis- 
sion, and very low excitation ionized gas tracers speak to 
the unusual environment of the region (Appleton et al. 
2006). This setting may be one in which the interstellar 
medium of young clusters is rapidly stripped, and thus 
young clusters do not have cocoons of gas to photoion- 
izc. However, though the young, compact star clusters 
are apparently not strong Ha emitters, optical spectra 
of bright H n regions do show prominent emission lines 
(e.g., Mendes de Oliveira et al. 2004). 

On the kiloparsec scales of complexes, HCG 31 hosts a 
population with masses that more closely resemble those 
found at higher rcdshifts. This suggests that, similar to 
galaxies at z = 1-2 with similar absolute magnitudes, the 
disks of galaxies in HCG 31 have higher turbulent veloci- 
ties than typically found in the local universe. This most 
plausibly arises because of the ongoing multigalaxy inter- 
actions in the system. In this sense therefore, the larger- 
scale star formation structures are reminiscent of those 
seen during the epoch of morphological galaxy transfor- 
mations in dense environments. 

4.2. Dwarf Galaxy Candidates 

Previous optical photometry suggested that galaxy F, 
composed of two separate H I peaks within the same 
cloud, is comprised of only a young population and that 
given its location in the tidal tail of the A+C complex it 
is likely a recently formed structure undergoing its first 
burst of star formation (Johnson & Conti 2000). Richer 
et al. (2003) find that galaxies E and F have similar oxy- 
gen abundances to A+C, verifying that they have been 
tidally removed from the A+C complex. The stellar 
populations are younger than the travel time from the 
A+C complex, and clearly are forming in situ (Mendes 
de Oliveira et al. 2006). 

Amram et al. (2004, 2007) found a flat / shallow rota- 
tion velocity in E and F using optical Fabry-Perot data, 
and thus classify them as tidal debris though dwarf galax- 
ies are expected to be slow rotators (Mendes de Oliveira 
et al. 2006). While E and H are merely fragments which 
are likely to fall back into A+C, F and R have H I column 
densities and Ha luminosities consistent with being true 
dwarf galaxies, and given their positions are less likely to 
fall back into the A+C complex. 

From the colors and luminosities of the SC candidates 
and complexes, the total stellar mass in F is > 10 6 M Q . 
While this is within the range found for dwarf galaxies, 
it is unclear whether it is sufficient to establish F as an 
independent entity within this compact group. From our 
analysis of the faint, diffuse emission in F, we find no 
evidence for an older stellar population, consistent with 
the Johnson & Conti (2000) analysis. 

4.3. Future Evolution 



The multiple interactions of low mass galaxies that 
characterize the HCG 31 system are evocative of hier- 
archical structure formation in process, and so we con- 
sider the likely end result, and compare it to early type 
systems in the local Universe that presumably under- 
went this type of dynamical processing at earlier times 
(z > 1). We thus examine if HCG 31 can be cast ac- 
curately as a local analogue of high-redshift hierarchical 
structure formation. 

At present, the sum of H I in HCG 31 as a whole is 
consistent with that expected from its member galaxies 
based on their magnitudes and morphologies (Williams 
et al. 1991; Verdes-Montenegro et al. 2001), though a 
large fraction of that gas has been pulled out from the 
individual galaxies into tidal structures. This suggests 
that HCG 31 has recently contracted as the gas has not 
yet been significantly reduced by a long period of star 
formation. The low radial velocity dispersions of the five 
primary galaxies (A, B, C, G, and Q) and the projected 
lengths of the H I tails are generally consistent with a 
bound system (Mendes de Oliveira et al. 2006). It seems 
quite likely that the individual galaxies in this group will 
ultimately merge into a single elliptical, as predicted by 
Rubin et al. (1990), in a wet merger process. Subse- 
quently, it is not likely to move along the red sequence 
to higher luminosities as there are no candidates for fu- 
ture dry mergers in the neighborhood. We then consider 
whether this system will become a so-called 'fossil group': 
the putative remnant of a coalesced compact group. 

A key component of the fossil group definition is a hot, 
extended X-ray medium with Lx > 0.5 x 10 42 /if 2 erg s _1 
(Jones et al. 2003). At present, HCG 31 has little ex- 
tended X-ray emission. Unlike in galaxy clusters and 
massive groups, the gravitational potential of HCG 31 
given its velocity dispersion is not sufficient to virialize 
gas to X-ray emitting temperatures. Therefore, the only 
possible source of heating is star formation whereby the 
mechanical energy from supernovae and stellar winds is 
thermalized into ~ 10 6 K gas. We can make a crude esti- 
mate of the potential X-ray luminosity (excluding com- 
pact sources) generated by conversion of the available 
reservoir of neutral gas into stars. In the following ar- 
gument, we use the numbers from Mihos & Hernquist 
(1996) for the efficiencies of gas-to-star and supernova- 
to-kinetic energy conversions. If 75% of the 1.6 x 10 10 M Q 
of H I is converted to stars, then ~ 9 x 10 7 supernovae are 
expected assuming the initial mass function of §2.5. If 
10~ 4 is the fraction of supernova energy converted to ki- 
netic energy, and all of that energy is ultimately radiated 
as X-rays, the integrated X-ray energy of all supernova 
(assuming E$n = 10 51 ergs) is 9 x 10 54 ergs. Over the 
course of a 150 Myr star formation episode, the power 
would thus be 2 x 10 39 erg s _1 , or approximately 0.4% 
of the minimum for a fossil group. Empirically, the dif- 
fuse emission from star formation observed locally is 1- 
5xl0 39 erg s^VO 1 © yr" 1 ) (Owen & Warwick 2009) , and 
so this estimate is not unreasonable. The total dynamical 
(including dark matter) mass of HCG 31 from H I veloci- 
ties (assuming spherical symmetry and dynamical relax- 
ation; Verdes-Montenegro et al. 2005) is 2 x 10 11 M , 
approximately two orders of magnitude smaller than the 
dark matter halo mass of the lower mass (and lower X- 
ray luminosity) fossil groups (e.g., Diaz-Gimenez et al. 
2008). Therefore, HCG 31 is a plausible candidate for a 
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low mass fossil group, where fossil group here applies to 
the dynamical history of a system rather than the spe- 
cific X-ray luminosity criteria as defined by Jones et al. 
(2003). 

We can also consider the current population of globular 
clusters, and investigate how that compares to those of 
low mass ellipitical galaxies in the local universe. While 
our value of 5n = 0.6 for a galaxy with My = —20.7 
(Ly = 1-6 x 10 10 L Q ) would be low for an elliptical of this 
absolute magnitude, late type galaxies have lower mass- 
to- light ratios, and therefore lower values of 5n even with 
the same number of clusters per unit stellar mass. A 
globular cluster system studied previously in NGC 6868, 
an E2 galaxy in the Tclcscopium group [a ~ 250 kms -1 ) 
with evidence of somewhat recent merger activity, has 
5n = 1.1, a value on the low end for ellipticals (Da Rocha 
et al. 2002a). Could such a galaxy be the result of a 
coalesced HCG 31? 

To compare globular cluster systems in galaxies of dif- 
ferent type, it is more convenient to use the T parameter, 
defined as the number of old globular clusters per unit 
stellar mass (Zepf & Ashman 1993). Assuming values 
of M/Lv=10 and 3 for ellipticals and late-type spirals, 
respectively, we obtain T = 1.7 ± 0.8 for HCG 31 as- 
suming a total globular cluster population of 111 ± 51 
for the entire system and a mass of 6 x 10 10 M . The 
corresponding T value for NGC 6868 (My = -21.9; 
Ly = 5 x 10 10 L ; M = 49 x 10 10 M Q ) using the globular 
cluster numbers from Da Rocha et al. (2002a) is 2.2± 1.1. 
Using My = -22.17 (Forbes et al. 2007) for NGC 6868 
gives T = 1.7 ± 0.9. Therefore, correcting for the differ- 
ing mass-to-light ratios of the two systems makes them 
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fully consistent within the uncertainties. Note that the 
T value for HCG 31 would likely be a lower limit, given 
that it is currently forming new star clusters, some with 
masses > 10 5 M Q (see Fig. 10); these massive clusters 
would only increase the T value if they survive the next 
few Gyrs intact. 

In sum, the likely outcome is that HCG 31 will have 
too low a mass and insufficient X-ray luminosity to meet 
the Jones et al. (2003) criteria for classification as a fossil 
group at z = 0. However, in terms of dynamical evolu- 
tion and Lx vs. halo mass, it could plausibly extend 
the relation seen in fossil groups to lower masses. The 
globular cluster populations of local ellipticals of simi- 
lar mass to HCG 31, however, are consistent with the 
current population of globular clusters in HCG 31 even 
if the ongoing star formation does not boost their num- 
bers. HCG 31 therefore does enable us to investigate hi- 
erarchical structure formation - the formation of a single 
elliptical from several smaller and gas-rich galaxies - in 
a local setting, though the likely result will be a smaller 
version of the end products of putative group coalescence 
that currently exist in the local Universe. 
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TABLE 1 

Star Cluster and Globular Cluster Candidate DATA a 



Region 


^B 435 b 


B - V° 


N d 


^SCCs 


/neb 


M v 606 (brightest) 8 


N G cc s h 


Sat 1 


A+C 


-20.0 


0.03 


11 


178 


0.29 


-15.5 


50 ±26 


0.5 ±0.3 


B 


-18.5 


0.17 


11 


105 


0.23 


-13.4 


39 ±22 


1.3 ±0.7 


E 


-16.5 


-0.03 


3 


17 


0.35 


-11.6 


(3) 


(0.3) 


F 


-15.8 


-0.08 


4 


16 


0.75 


-11.4 


(3) 


(0.5) 


G 


-18.9 


-0.01 


11 


86 


0.33 


-13.4 


20 ± 13 


0.6 ±0.4 


H 









3 


0.67 


-11.2 









a The magnitudes in this table have been corrected for Galactic extinction only. 
b Absolutc magnitude for the entire region. 

c The B — V color for each region is reproduced from Table 2 of Lopez-Sanchez et al. (2004). 

d The number of complexes in each region (see §2.3). 

e The number of star cluster candidates in each region (see §2.1). 

'The fraction of star cluster candidates in each region with colors consistent with nebular 
emission (see §2.1 and the lefthand panel of Fig. 3). 

g The absolute magnitude, My 606 , of the brightest SC candidate in each region. The mag- 
nitude has not been corrected for intrinsic extinction. 

h The total calculated number of globular cluster candidates after correction for incomplete- 
ness (see §3.4). The numbers in parentheses were estimated from only one object. 

'The specific frequency of globular cluster candidates (see §3.4). The numbers in parentheses 
were calculated from only one object. 



TABLE 2 

Star-forming Complexes 



Region 


Complex 


RA 




Dec 


M V eoe B 435- 


-^606 


V606 _ ^814 


Diameter (pc) 


A+C 


1 


5:01:39.7676 


-4: 


15:11.865 


-14.09 


0.31 


0.23 


539 


A+C 


2 


5:01:38.6204 


-4: 


15:14.205 


-13.25 


0.25 


0.47 


331 


A+C 


3 


5:01:38.4399 


-4: 


15:13.545 


-13.88 


0.26 


0.16 


452 


A+C 


4 


5:01:38.0909 


-4: 


15:24.225 


-14.75 


0.20 


0.26 


414 


A+C 


5 


5:01:38.7127 


-4: 


15:32.985 


-15.37 


0.47 


0.64 


306 


A+C 


6 


5:01:38.2434 


-4: 


15:33.285 


-15.20 


0.41 


0.38 


271 


A+C 


7 


5:01:37.7319 


-4: 


15:32.715 


-15.81 


0.85 


-0.38 


198 


A+C 


8 


5:01:37.6858 


-4: 


15:32.745 


-14.71 


0.60 


-0.09 


191 


A+C 


9 


5:01:37.8062 


-4: 


15:27.945 


-17.73 


0.51 


-0.21 


489 


A+C 


10 


5:01:37.7500 


-4: 


15:28.965 


-15.77 


0.36 


-0.00 


265 


A+C 


11 


5:01:37.7019 


-4: 


15:31.305 


-17.01 


0.27 


-0.04 


418 


B 


1 


5:01:34.7918 


-4: 


16:00.254 


-14.34 


0.38 


-0.12 


288 


B 


2 


5:01:35.3112 


-4: 


15:47.414 


-13.71 


0.32 


0.39 


274 


B 


3 


5:01:35.2350 


-4: 


15:50.534 


-13.06 


0.47 


0.20 


166 


B 


4 


5:01:35.5840 


-4: 


15:41.834 


-13.52 


0.51 


0.18 


265 


B 


5 


5:01:35.8327 


-4: 


15:40.334 


-13.84 


0.62 


-0.21 


265 


B 


6 


5:01:35.7685 


-4: 


15:40.154 


-13.00 


0.43 


0.10 


207 


B 


7 


5:01:35.8327 


-4: 


15:42.614 


-12.54 


0.77 


-0.31 


133 


B 


8 


5:01:35.9169 


-4: 


15:42.314 


-13.34 


0.31 


0.18 


257 


B 


9 


5:01:35.7966 


-4: 


:15:44.114 


-12.50 


0.27 


0.42 


166 


B 


10 


5:01:34.9863 


-4: 


15:55.874 


-13.08 


0.34 


0.47 


248 


B 


11 


5:01:35.2150 


-4: 


15:51.974 


-13.57 


0.35 


0.34 


208 


E 


1 


5:01:37.7058 


-4: 


15:57.405 


-13.47 


0.40 


0.04 


299 


E 


2 


5:01:37.5314 


-4: 


15:56.715 


-15.80 


0.20 


0.36 


875 


E 


3 


5:01:37.2887 


-4: 


15:48.075 


-12.89 


0.36 


-0.00 


277 


F 


4 


5:01:39.8016 


-4: 


16:21.525 


-15.11 


0.83 


-0.63 


550 


F 


5 


5:01:39.8598 


-4: 


16:22.515 


-13.60 


0.38 


-0.24 


333 


F 


6 


5:01:40.1947 


-4: 


16:26.835 


-14.27 


0.67 


-0.47 


373 


F 


7 


5:01:40.2589 


-4: 


16:28.755 


-12.48 


0.36 


-0.30 


286 


G 


1 


5:01:44.0455 


-4: 


17:21.045 


-14.68 


0.31 


0.41 


252 


G 


2 


5:01:43.9492 


-4: 


17:20.985 


-14.83 


0.34 


0.41 


298 


G 


3 


5:01:43.9212 


-4: 


17:22.665 


-13.89 


0.45 


-0.01 


190 


G 


4 


5:01:43.8610 


-4: 


17:22.305 


-13.86 


0.37 


-0.01 


187 


G 


5 


5:01:43.7687 


-4: 


17:19.245 


-13.98 


0.36 


0.21 


214 


G 


6 


5:01:43.7206 


-4: 


17:19.485 


-13.06 


0.17 


0.15 


139 


G 


7 


5:01:43.6685 


-4: 


17:17.625 


-15.05 


0.63 


-0.23 


307 


G 


8 


5:01:43.8289 


-4: 


17:14.445 


-15.30 


0.38 


0.17 


407 


G 


9 


5:01:43.9011 


-4: 


17:14.085 


-14.99 


0.49 


-0.00 


357 


G 


10 


5:01:43.9492 


-4: 


17:15.105 


-13.50 


0.32 


0.37 


198 


G 


11 


5:01:44.0816 


-4: 


17:13.965 


-13.58 


0.26 


0.20 


184 
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TABLE 3 

Ultraviolet Fluxes, Infrared Luminosities, and Star Formation Rates 







Swift 




GALEX 


Spitzer 




UV+IR 


U 


UVM2 UVW1 


UVW2 


NUV 


FUV 




SFR 


Region 




(10- 14 erg s- 1 


cm' 2 A- 1 ) 






(log(crg s- 1 Hz" 1 )) 


(M© yr- 1 ) 


A+C+E 


2.1 ± 0.1 


4.8 ±0.2 3.0 ±0.1 


5.2 ±0.2 


4.3 ±0.4 


7.9 ±0.8 


30.29 


8.11 


B 


0.34 ± 0.01 


0.81 ± 0.03 0.46 ± 0.02 


0.89 ±0.04 


0.77 ±0.07 


1.3 ±0.1 


28.84 


0.78 


F 


0.25 ± 0.01 


0.31 ±0.01 0.31 ±0.01 


0.27 ±0.02 


0.16 ±0.02 


0.29 ± 0.03 


28.06 


0.19 


G 


0.56 ± 0.02 


1.58 ±0.05 0.82 ±0.03 


1.74 ±0.08 


1.5 ±0.1 


3.0 ±0.3 


29.21 


1.47 


Q 


0.033 ± 0.001 


0.079 ± 0.004 0.036 ± 0.002 


0.088 ± 0.008 


0.098 ± 0.009 


0.17 ±0.02 


27.46 


0.07 


Total 












30.34 


10.62 
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Fig. 1.— A five-color composite image of HCG 31 composed of purple (NUV GALEX), blue (B 435 HST), yellow-green (V 60 6 HST), 
orange (8 fim Spitzer), and red (/gi4 HST) images; see § 2.8. Group galaxies and tidal debris are labeled in light gray; HCG 31D (the 
red spiral in the upper right) is in the background. Young, massive stars are evident in blue and purple, and dust emission shows up as 
light orange. In general, the brightest nebular emission regions glow pink, because they have both strong UV and dust emission. However, 
regions such as tidal feature F with less pronounced IR emission appear turquoise due to H a emission falling within the wide Veoe- The 
30" scale corresponds to roughly 8 kpc. North is up, and east is to the left. 




Fig. 2. — ACS three-color image of HCG 31. The gray polygons denote the boundaries of the physical regions associated with the galaxies 
and tidal structures in HCG 31 as labeled. They were defined by eye to encompass the optical light associated with each region. 
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Fig. 3. — B435-V6O6 vs. Veoe~ -^814 color-color plots of the star cluster (SC) candidates (black crosses) with Mv 606 < —9 and complexes 
(open purple diamonds) for the entire ACS field of view of HCG 31. The large black cross in the lower right corner indicates the median 
errors of the SC candidates; a Galactic reddening vector of ^4606=1 is overplotted in the upper right corner. The solid, red curve is a 
Bruzual & Chariot (2003) evolutionary track for an instantaneous burst with 1/5 solar metallicity and a Salpeter initial mass function. 
The dashed green curve is an evolutionary track which includes nebular emission as described in section §2.5; the large numbers of SC 
candidates and complexes to the left of this track have colors indicating ages < 5 Myr. The positions of many of them below the dashed 
line is consistent with internal reddening shifting their colors. From visual inspection, most of the identified SC candidate and complexes 
have colors that are consistent with 10 7 -10 8 yrs. For these sources, the tracks are almost parallel to the reddening vector, and there is 
considerable age-reddening degeneracy in this region of color space. Few globular cluster candidates are expected in this figure given the 
bright V606 luminosity limit. 



AC 





Fig. 4. — Star-forming complexes are shown in the three-color ACS images of HCG 31. The whole field is shown in the lower left, with 
the galaxies and tidal structures labeled. Close-ups of A+C and F are shown in subsequent frames. Measured complexes are indicated by 
solid circles with colors chosen to enhance visibility. Note that these complexes may encompass several star clusters. The thick, gray bar 
in each image represents 1 kpc. 
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Fig. 5. — Same as Figure 4 for B, E, and G. The thick, gray bar in each image represents 1 kpc. 
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Fig. 6. — A V606 vs. V606 - ^814 color-magnitude diagram for the SC candidates (black crosses) with Mv 608 < —9 and complexes (open 
purple diamonds). The evolutionary tracks for 10 5 , 10 6 , and 10 7 Mq clusters are indicated with the curves as labeled; dashed curves 
include nebular emission. Galactic extinction of ^4606= 1 has been indicated with a labeled vector, and log(ages) in years are marked. The 
median color-magnitude error is indicated by the large black cross in the lower right corner. SC candidates consistent with having masses 
> 1O B M0 are found; the majority have V606 - ^814 colors ranging from 0.2-0.6. 
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Fig. 7. — In each panel, the Mv 606 distributions of star cluster (SC) candidates with My 606 < —9 for each region as labeled. Filled black 
histograms represent the distribution for SC candidates with nebular colors (see caption to Figure 3) . The histograms overall are consistent 
with power-law distributions as found in other star-forming systems. 




Fig. 8. — The absolute magnitude, Mv 606 , of the brightest cluster in each region vs. log(TV) (where N is the number of star cluster 
candidates brighter than Mv 606 = —9 which excludes the majority of globular cluster candidates) for the entire ACS field. Each region 
is as labeled in the legend. Though no color criteria were applied to the selection of the brightest clusters, they are all young and blue. 
The errors in the ordinate are from poisson counting statistics, and the abscissa errors arc photometric. The SC candidate populations in 
HCG 31 are consistent with the best-fitting linear fit (solid line) from the study of Whitmorc et al. (2007) of 40 local star-forming galaxies, 
none of which were in compact groups. The scatter is also consistent; the dashed lines encompass 90% of their sample. Whitmore et al. 
(2007) interpret this correlation as a consequence of larger SC populations being statistically more likely to sample the bright end of the 
luminosity function. 
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Fig. 9. — B435— V606 vs - V606 _ ^814 color-color plots of the SC candidates with My eae < —9 and complexes in HCG 31. In each panel, 
black dots indicate all SC candidates from the entire ACS field of view, black crosses represent candidates within the region as labeled 
and shown in Figure 2. Larger, blue crosses are the more luminous (My eoe < —11) SC candidates, the most luminous SC candidate is 
indicated with an open, green star, and open purple diamonds are the complexes. The large black cross in the lower right corner indicates 
the median color errors for the SC candidates in each region. Evolutionary tracks and reddening vectors are as described in the caption 
to Figure 3. The SC clusters are concentrated in the galaxies A+C, B, and G, though some are also found in the tidal features E, F, and 
H. Complexes are found in both galaxies and tidal features, with the exception of H. From comparison of the colors with the evolutionary 
tracks, all regions have SC candidates with ages < 10 7 yr and significant nebular emission, and only F and H show no evidence for SC 
candidates with ages > 500 Myr. 
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Fig. 10. — V606 vs - Vfeoo - ^814 color-magnitude plots of the SC candidates with MV606' 
black dots indicate all SC candidates from the entire ACS field of view, black crosses represent candidates within the region as labeled and 
shown in Figure 2, the open black star is teh most luminous SC candidate, and open purple diamonds are the complexes. The median 
color-magnitude error in each region is shown with a large black cross in the lower right corner. Evolutionary tracks and reddening vectors 
are as described in the caption to Figure 6. By comparison with the evolutionary tracks, A+C has both complexes and SC candidates that 
are consistent with masses up to 10 Mq, notably more massive than those in G and B. The tidal features E and F also host massive 
complexes, though their SC candidate populations are not as massive as those found in the galaxies A+C, B and G. 
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Fig. 11. — Color-color diagram for all point sources with —13.8 < My 606 < — 7.8 (20 <V606< 26) where the strict point source selection 
criteria are defined in §2.1; this selection is designed to identify globular cluster candidates. The red dotted region denotes the expected 
colors for Milky Way globular clusters based on the Harris (1996) catalog; see §3.4 for details. The solid, red line within the detection region 
shows the predicted colors for old stellar populations from the simple stellar population models of Maraston (2005). Black symbols indicate 
the regions where the point sources are located (see legend for key); green circles underlie the 43 objects whose photometric error ellipses 
intersected the GC selection region. The median color-color errors are plotted as a large cross in the lower left corner. For reference, the 
nebular (blue, dashed) and normal (blue, solid) evolutionary tracks as described in §2.5 and an Aqqq=1 reddening vector are overplotted. 
Note the lack of strictly defined point sources with nebular colors compared to Figure 3; this shows that the vast majority of nebular 
sources are marginally resolved in the V606 image. 
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Fig. 12. — Raw V606 luminosity function (histogram) for all GC candidates with no correction made for any (small) background 
contamination or for photometric completeness. The dotted line shows the expected luminosity function for GCs located at the distance 
of HCG 31: a peak magnitude of My = -7.4 (corresponding to Vaoa= 26.3) and a dispersion of a = 1.2. The plotted function is not a 
fit to the histogram, but is shown to indicate its expected shape. 
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Fig. 13. — A Vega-magnitude calibrated, point-source subtracted, boxcar-smoothed (with a 13 pix window), B435-/814 image of HCG 31. 
The image has been corrected for Galactic extinction. The scale is linear, and plotted colors range from —0.9 (black) to 1.9 (white). 
The blue-red transition occurs at -B435— 7gi4=0.4, and -B435— 7gi4=0.7 and 1.0 contours are plotted in green. Regions are labeled as in 
Figure 2. The lack of red in much of the F region indicates that there is no evidence for an older stellar population in the diffuse emission 
in that region. In contrast, the tidal material pulled out to the north of A+C and coincident with E and H clearly has redder colors 
(5435-^814=0.4-0.9). 
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Fig. 14. — V606 vs - B435— Jgi4 color-magnitude diagram of the SC candidates and complexes; symbols and tracks are as described in the 
caption to Figure 6 with the addition of large crosses and diamonds to mark the objects in F. The vertical lines at -B435— /si4=0.4, 0.7, and 
1.0 mark the blue-red transition and the green contour levels in the image shown in Figure 13. S435— Jsi4 color is not sensitive to nebular 
emission, and old (> 500 Myr) star clusters are expected to have S435— i8i4> 1.0. From their colors, the SC candidates and complexes in 
F are notably young with B435— Jsi4 values extending from —0.3 to 0.7. 



